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ABSTRACT 

Polars (or AM Her systems) are cataclysmic variables without a disc due to the strong 
magnetic field of the white dwarf. Most of their emission comes from the region where the 
accretion column impacts the white dwarf and cools through cyclotron and bremsstrahlung 
processes. We present a new code, CYCLOPS, to model the optical emission from these 
systems including the four Stokes parameters. It considers a three-dimensional region with 
the electronic density and temperature varying following a shock-like profile and a dipolar 
magnetic field. The radiative transfer is solved in steps considering the solution with non- 
null input radiation. The footprint of the column in the white-dwarf surface is determined 
by the threading region in the equatorial plane, i.e., the region from where the flow follows 
the magnetic lines. The extinction caused by Thomson scattering above the emitting region 
is optionally included. The fittings of observational data are done using a hybrid approach: a 
genetic algorithm is used to seek for the regions having the best models and then an amoeba 
code refines the search. An example of application to multi-wavelength data of V834 Cen 
is presented. The fit found is consistent with previous parameters estimates and is able to 
reproduce the features of V834 data in three wavebands. 

Key words: novae, cataclysmic variables - radiative transfer - methods: numerical - polari- 
sation - magnetic fields 



1 INTRODUCTION 

Cataclysmic variables (CVs) are short-period binaries consisting of a white dwarf (primary) and a late-type main-sequence star. The sec- 
ondary star fills its Roche lobe, losing material to the primary by the inner Lagrangian point, L\. Polars, also called AM Her systems after 
their prototype, are CVs in which the magnetic field reaches the 10 to 200 MG range in the primary surface. Such a strong field significantly 
affects the dynamics of the mass transfer preventing the formation of the accretion disc, a typical feature of non-magnetic CVs. In polars, 
the material from L\ follows a ballistic trajectory on the orbital plane (horizontal stream) up to the magnetic coupling region, which is also 
called the threading region. From this region onwards, the gas-flow effectively traces the geometry of the magnetic field, forming an accre- 
tion stream which reaches the white dwarf (WD) with supersonic velocity and ionised. Another consequence o f the strong mag neti c field is 



the sy nchronization of the white-dwarf rotation with the orbital revolution. Reviews on polars can be found in lCroppen d 19901) and lWamei 



lis 



(19951). 

Near the white dwarf, the magnetically controlled flow creates a shocked region that emits most of the system flux in optical and X-ray 
wavelengths. It is usually called the post-shock region. Its location on the white-dwarf surface depends on the geometric properties of the 
magnetic field and the threading region distance. Considering a dipolar magnetic field, the base of the accretion stream could not be at the 
magnetic pole since the threading region would be too far. Specifically, given a magnetic field geometry, the distance from the threading 
region to the white dwarf defines the angle between the magnetic pole and the post-shock region (Cropper'' 1989'). 

The optical flux and polarisation of polars are highly modulated with the orbital phase. These observations can be explained by the 
cooling of the post-shock region through cyclotron process. This emission depends strongly on the viewing angle and the physical properties 
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of the emitting region. Stokes tecliniques are consequently a powerful tool to diagnose the physical and geometrical properties of the accreting 
region. 



An important first step in the modelling of polar observation s was given in '80s by different groups jChanmugam & Dulklll981 : 
Mepgitt & Wickramasingha 1982 : Wickramasinghe & MeggitJl 19851) . These works assumed the emitting region as a point source with fixed 



values of magnetic field induction, electronic temp erature, and size parameter, A - this last parameter carries information about the electronic 
density and the optical depth of the emitting region. IWickramasinghe & Meggittl ( ll985h present a grid of fluxes, and linear and circular polar- 
isations as function of viewing angle and harmonic for some models with different values of A and electronic temperature, which can be used 
to fit optical data of AM Her systems. However, there are many pieces of observational evidence of an extended emitting region in polars. A 
natural impr ovement therefore is a t wo-dimensi onal modelling using a sum of po int sources with distinct parameters to represent the emitting 
region (e.g.. iFerrario & Wickramasingheiil99ol : |Potter. Hakala & Croppei 1998h . This approach is based on an optically thin hypothesis in 
which the overlap of the emitting points do not produce moderate opac ities. The density and temperature can also vary in the radial d i rectio n 
as a consequence of th e shock. This was introduced in the models by Wickramasinghe & Ferrario 1 1988h . Wu & Wickramasinghe ( 199C ) 



and 
and 



Canalle et al 



Potter et al] 1 2004|). Calculations on the shock structure of magnetically-channelled accretion flows can be found in Cropper et al. 1 1999b 



12005). 

In this work, we present a code to simulate the continuum optical emission from polars. We call it CYCLOPS, an acronym for CY- 
CLOtron emission of PolarS. It incorporates the radiative transfer of the cyclotronic process in the post-shock region adopting a three- 
dimensional treatment and specific functions to describe the variation of physical quantities. For the first time, the radiative transfer is solved 
in a stepwise manner considering in each step the homogeneous solution of the transfer equation with non-null input radiation. The shape 
of the emitting region is defined by following the magnetic field lines from the threading region to the white-dwarf surface. A procedure to 
fit observational data is also included. In Sec.|2l we describe the model. Some examples of results and comparisons with previous works are 
shown in Sec. [5] An application to V834 Centauri's data is presented in Sec.E] The last section presents a summary of this work. 



2 A NEW MODEL TO THE EMISSION FROM POLARS 

Our objective is to model the cyclotronic emission from the post-shock region in the accreting columns of polars. It is considered as an 
extended three-dimensional (3D) region composed by thermal electrons immersed in a dipolar magnetic field. The code considers electron 
temperature and density profiles in the radial and tangential direction of the emitting region, i.e., they are not constant inside the emitting 
region. In this section, we describe our model as well as some details of its implementation, which has been done using as framework the 
Interactive Data Language (IDL). 



2.1 Defining tiie emitting region 

In this section, we present an overview of the adopted scenario to model the optical emission coming from polars, which is schematically 
represented in Figure[T] The geometry of the emitting region is defined considering a poloidal magnetic field that captures the material from 
the ballistic trajectory. This coupling region is called the threading region. Our model assumes that it is a 2D surface on the equatorial plane 
of the binary system. Below we enumerate the main definitions and assumptions of our model. 

Coordinate system. The system as a whole is defined by a spherical coordinate system centred in the white dwarf (WD) being the 
polar axis parallel to the rotation axis of the white dwarf (this is also the direction perpendicular to the orbital plane). A reference point 
inside the emitting region defines the longitude zero. We call it the anchor point, which is represented as a black dot in Figure [T] (right). The 
colatitude of the anchor point is (5. 

Inclination. The system can be seen by any inclination, which is defined by the angle between the WD rotation axis and the observer. 

White dwarf. The primary is described by an opaque sphere of arbitrary radius, Rwd- The eclipse of the emitting region by the WD 
is consequently included. From the definition of the coordinate system, the WD angular momentum has the same direction and sense of the 
polar axis. The default sense rotation of the WD occurs in the counter-clockwise sense as seen by an observer at a positive latitude. 

Magnetic field. The direction and intensity of the magnetic field in the emitting region follow the prescription of a centred dipolar 
field. The magnetic axis can have any angle (colatitude) to the rotation axis. Also its longitude is free. Thus, the rotation axis, magnetic axis, 
and accreting region are not necessarily in a same plane: this allows the magnetic field lines to depart considerably from the radial direction. 

Geometry of the emitting region. Its shape is of a magnetic tube whose inferior and superior ends are sections of spherical surfaces 
with radius [bRwo] and [{b + h)RwD], respectively. This means that a region with b = 1 and h — 0.1 runs from the WD surface to 10% 
Rw D in height. The borders are defined by the magnetic field lines that delineate the threading region in the equatorial plane of the binary 
system as shown in Figure [T] This procedure provides a better approach than an ad-hoc choice and a tool to, in future works, understand 
the threading region in polars, i.e., the region where the matter coming from the secondary is coupled to the magnetic field lines. This also 
prevents the choice of an emitting region resulting from an unlikely threading region. In order to calculate the "walls" of the emitting region, 
we first follow the magnetic field line that passes through the anchor point up to the crossing with the equatorial plane of the binary system. 
This defines a point so forth called the threading point at a radius of Rth and a longitude of longth'- it is represented as an asterisk in Figure 
[T](left). Then we consider a 2D threading region extending from (1 — Ar) Rth to (1 + Aj?) Rth in radius and 2Aiong in longitude. The 
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magnetic field lines that define this boundary are followed back to the WD surface, and so define the border of the emitting region. An 
example of the emitting region on the white dwarf is presented in Figure [T] (right). In this procedure, Rtu and longth are determined by 
P, the colatitude of the anchor point, and the magnetic field axis latitude and longitude. The extension of the threading region accounts for 
the dispersion in velocities and angular momentum of the accreting matter. We would like to note that the threading point is not necessarily 
at the centre of threading region, which extends from \longth — 2/; A;o„g] to [longth + 2(1 — /i)A;o„g] in longitude. In this manner, the 
threading point location is controlled by the parameter /; and can be at any longitude between its least and maximum values. This brings 
about an anchor point that can be off centre too. 

Electronic density and temperature. Our code allows us to adopt any analytical function to represent the matter density and 
temperature in the emitting region. In this work, we have considered the following options for their variation. In the tangential direction, no 
variation at all or an exponential decrease from the anchor point to the border of the region. As the anchor point is not necessarily centred, as 
explained above, the point of highest density/temperature can be displaced from the geometrical centre of the region. In the radial direction, 
the density decreases exponentially and the temperature increases exponentiall y from the WD surf ace to above. The used functions try to 
simulate the temperature and density shock solutions as shown in Figure 1 from lCropper et alj ( 119991) including cyclotron cooling. However, 
they can be scaled in value and in height. 

Phase. The zero phase corresponds to the anchor point in the meridian of zero longitude. Therefore, in this phase the anchor point is 
located in the plane containing the white-dwarf rotation axis and the observer direction. 



2.2 The radiative transfer 



In the previous section, we have described the main assumptions of our code relative to description of the emitting region. Now we turn to 
the details of the radiative transfer solution for the cyclotron emission. 

The radiative transfer is solved in a discrete approach as explained here. We initially define a Cartesian coordinate system in which the 
z axis is the observer direction, i.e., the direction of the radiation propagation. The x axis is in the orbital plane, and the y-axis completes 
the coordinate system. It is called the observer reference frame. A cube containing the emitting region is then defined and divided in N"^ 
voxels (elements of volume). In this representation, the emitting region is composed by (a maximum of) lines of sight. In each of them, 
the radiative transfer is solved in steps of homogeneous voxels. Specifically, the density, temperature, and magnetic field are allowed to vary 
from voxel to voxel, but are constant inside a voxel. The values assumed for the calculations are that of the central point of the voxel. The 
number of voxels, defining the resolution of the grid, are freely chosen to properly account for the variation of the parameters and geometry 
of the region. The Stokes parameters that represent the integrated emission from the region are the sum of Stokes parameters of all line of 
sights. The voxel structure is redone for each phase in which we run the code. 



The solution of the transfer of polarised ra diation in magneto-active plasma is found in lPacholczvk & Swiharll ( Il975l . eqs. 7 to 13) or 



Zheleznvakov. Suvorov & Shaposhnikovl dl974l) . It considers a non-null incident radiation and therefore should be used here, because from 



the second voxel and forth the emission from the preceding voxel should be considered as input. This gives us confidence that we have 
an adequate treatment to any optical depth, even those intermedia t e from the optically thin to thick regimes. In Appendix [A] we present 
an alternative expression to the solution of Pacholczvk & Swihart i 19751) and some misprint s in the original formu l ae. In order to avoid 
computational problems in voxels with optical depth much smaller than 1, the formalism of IPacholczvk & Swiharj ( Il975h for negligible 
absorptio n was used (their eg. 14 ) . 

The Pacholczvk & Swiharll 1 1975)'s solutions are defined in a coordinate system having one axis parallel to the component of the 
magnetic field perpendicular to the radiation propagation direction. We call it the plasma reference system. The magnetic field direction 
changes from voxel to voxel and thus the plasma reference sys tem. Consequently this system is not appropriated to represent the emergent 
emission from the region as whole - see discussion presented in lPacholczvk & Swiharll ( 1 19671) . The incident and emergent Stokes' vectors of 
each voxel are therefore represented in the coordinate system of the observer defined above. As a consequence, to solve the radiative transfer, 
these vectors have to be transformed to and from the plasma reference system in each voxel . 

The cyclotronic emission coefficients have been calculated using t he expressions from iSeard & Bakej ( il962h and ichanmuganij \\9%^ . 
Our results agree very well with Figure 1 from 'chanmugam (1980). The absorption c o efficie nts have been obtained from the Kirch- 
hoff's Law, consi dering the thermal regime, a s proposed by Meg gitt & Wickramasinghd ( Il982r) . The particular solution of the radiative 
transfer equation I Pacholc zvk & Swihart 19751) in this case equals to the Planck function (see Equation IA4t . Following the procedure of 
Iwickramasinghe &^eggi'tl' /l985). the b remsstrahlune abso rption has also been considered. The Faraday mixing coefficients can be found 
in IMeggitt & Wickramasinghe (.1982) or IPacholczvkl (19771). In this work, we have not included any scattering process in the radiative 
solution. 

The computation is done considering a fixed frequency. Consequently, the harmonic number of the cyclotron emission varies throughout 
the region, because the magnetic field changes in space due to the dipolar geometry. When fitting observational data, it must be taken into 
account the width of the observing bandpass. We will return to this topic in Section|4| 

In order to test the algorithm of integration by steps of the radiative transfer, we constructed regions with homogeneous physical proper- 
ties (magnetic field intensity and direction, electronic density, and tempe rature). The resulting spectra of the ste pwise solution, for different 
choice s of parameters, was confronted to the one-step solution provided in Meggitt & Wickramasinghe 1 1982 ) and Wickramasinghe & Meggit3 
19851) . The two solutions are in very good agreement. 
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2.3 Extinction by the upper portion of the accreting column 

The magnetically controlled flow in polars begins at the threading region and ends at the WD surface. Our current view of polars assumes 
that only a portion with enhanced density near the WD, the post-shock region, emits by the cyclotronic process. There is, however, material 
in the flow outside the emitting region that could, in principle, attenuate this emission, if it is in between the observer and the post-shock 
region. This portion of the accretion column will be called the attenuation region. We have, then, included in our code the calculation of this 
region for which an extinction can be assumed. The attenuation region is divided in voxels and limited by the same magnetic field lines that 
delineate the emitting region, thus a given line of sight (from the 7V^ ones that represent the emitting region) has above a specific number 
of voxels to which some attenuation can be assigned. All the voxels have the same value of attenuation, hence each line of sight has a total 
attenuation proportional to the number of voxels above it. The inclusion of this attenuation is optional in our code. If the user decides to use 
it, the program chooses the best value that fits the observed data (see next section). Each voxel has a maximum attenuation limited by the 
Thomson scattering cross section and the minimum density in the emitting region, which is achieved in its superior border. 

2.4 Fitting observational data 

Understanding the optical emission of AM Her systems can be divided in two main tasks: (1) to propose a model to the system and to calculate 
its emission; (2) to decide, among the many possible sets of model parameters, the one that best fits the observations. This section concerns 
the procedure to fit the data. Our code aims to reproduce the total, linear and circular polarised fluxes, and the angle of the linear polarisation 
as function of the orbital phase. Our code produces images of the four Stokes parameters (see Section O. However, the observational data 
have not enough angular resolution to see the post-shock regions as extended objects. So the code integrates the images in one speficic phase 
providing a point per Stokes parameter. 

A model is described by a set of parameters which are explained in detail in Section[2T| The CYCLOPS parameters are presented as a 
quick list below. 

• i, inclination; 

• /3, the colatitude of the emitting region; 

• b, the distance from the WD centre to the base of the region in units if WD radius. In spite of be a free parameter in the code, we assume 
it equal to one in all calculations presented in this work; 

• h, the length of the emitting region in the WD radial direction; 

• AT?, one-half of the fractional radial extension of the threading region; 

• l^long, one-half of the azimuth extension of the threading region; 

• /i, this parameter controls the position in the longitudinal direction of the threading point relative to the geometrical centre of the 
threading region. If the value of /; is (0.5; 1.0), the threading point is at the border of smallest longitude (centre; border of largest 
longitude) of the threading region. We would like to recall that this parameter controls the point of maximum density and/or temperature in 
the emitting region; 

• B, the magnetic field intensity in the magnetic pole on the white-dwarf surface from which we calculate the magnetic moment of the 
dipole; 

• Biat, the latitude of the magnetic axis; 

• Biong, the longitude of the magnetic axis; 

• Ne, the maximum value of the electronic density; 

• T, the maximum value of the electronic temperature. 

The user must specify the frequency(ies) in which the data have been collected and is(are) not a free parameter. 

There are parameters that are not part of the physical description of the system, but should be specified in the calculation. They are the 
number of orbital phases and the number of voxels in one dimension in the cube inscribing the emitting region (see Section lZTt . 

When fitting an observational data set, there are five quantities that are calculated for each model. We discuss each one in the following 
paragraphs. 

SCALE FACTOR TO THE CYCLOTRON FLUX. The code firstly calculates the scale factor that multiplies the circular polarised flux of the 
model in order to provide the best match with the observed data. This flux component can only be explained by cyclotron emission and is, 
therefore, the best choice to fix the scale factor of the model. The same factor is then applied to the total and the linear polarised flux of the 
cyclotron model. 

THE CONTRIBUTION OF AN ADDITIONAL UNPOLARISED SOURCE OF LIGHT. After fixing the scale of the model flux using the circular 
polarisation data, the code calculates the flux of an unpolarised and constant source of light that should be added to the total flux of the model 
to best fit the observed flux. The addition of such a source of light is usually necessary when fitting data and can be easily justified by the 
contribution of other components of the system, as the white dwarf, for instance. 

AN OFFSET IN PHASE. The code also calculates the offset between the observed and model phases that best reproduces the data. 

AN OFFSET IN POLARISATION ANGLE. When fitting the linear polarization, we included an offset in polarisation angle to obtain the best 
agreement between model and data. 
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THE MAXIMUM EXTINCTION IN THE ATTENUATION REGION. As explained in Section [23] our code can include the extinction caused by 
the upper portion of accreting column in front of the emission from the post-shock region. If this option is set, the attenuation is fitted with 
the constraint of the maximum value it can assume for Thomson scattering in the attenuation region. 



To search for the post-shock region properties in the domains of such a la rge number of input parameters , we have adopted a two-step 
approach. Initially, we use the pikaia implementation of a genetic algorithm bv lCharbormeaul ( 1199511^1 to look for the region in the parameter 



space associated with the maximum of a figure of merit related to the goodness of the fitting (g.o.f.). We then use an amoeba code (Press et alj 
1992h to refine the search for the best model. The g.o.f. is measured by the or a modified expressions similar to it. We have implemented 
to weight or not each parcel of the with the observational error The g.o.f. of a model is the sum of the for each observed curve: the 
total flux, the linear polarised flux, the circular polarised flux, and the position angle of the linear polarisation. For two or more bandpass 
data, the g.o.f. is obtained by the sum of the for each band. 



3 GENERAL RESULTS 

An important aspect of the emission from polars is the high variability with phase. The reason is two-fold. Firstly, the cyclotron radiative 
coefficient is strongly dependent on the viewing angle of the magnetic field. Second, the physical depth of the emitting region changes with 
phase in a region not spherically symmetric. A specific phase has, therefore, a total emission different from another phase as a result from 
matter distribution and view angle of the magnetic field. Both aspects of the cyclotron emission in polars is properly treated by our model. In 
this section, we show some results that illustrate the above discussion. 

CYCLOPS has as parameters the values of ele ctronic density and temperature. Previous models tend to use the size parameter, A. Its 
definition is fe.g.. lMeggitt & Wickramasing h3 ll982l) : 



A = 2.01 10^ 



V 105cm/ VlOi^cm^V 



iV, \ /3 10^G 



B 



As we are using a 3D treatment, a value of A does not specify a unique model: we have a family of models with different values of electronic 
density and height of the emitting region, h, that can produce a same value of A, supposing the magnetic field a s fixed. 



We would like to start the presentation of our results with their validation in comparison with those of iFerrario & Wickramasinghe 



who presented 2D extended models of polar s. Figure |2] (left) presents our model calculated with the parameters equivalent to the 



Figs. 2-5(c) from ^ Ferrario & Wick ramasinghe ('l990'). We have adopted constant values of electronic density and temperature without the 
addition of a component of constant and unpolarised flux. The curves represent the results for two values of h, 0.03 and 0.10 (A^e equ als to 
cm"'', respectively). These numbers provide the same A of the models of Ferrario & Wickramasinghe ( 1990l) for a 



radial line of sight passing in the centre of the emitting region. We have adopted a frequency of 6.6 lO'^ ^ Hz (optical wavelengths) in order 
to have the harmonic number near 8. Our curves are almost identical those fromj^errario & Wickramasinghe. (.199Q) . The small difference 
are acceptable, giving the different methodologies. The maximum optical depth (considering all the phases and all the line of sights) occurs 
around phase 0.5 and is 0.73 (0.23) for the model using h = 0.03 (0.10) and is exactly where our two models present distinct values (Figure 
(2] right). That shows that the differences quickly appear as the opacity enlarges, as expected, and the optical depth depends on phase. 

On the other hand, the differences between our approach and 2D calculations is illustrated in Figure [3] It presents the results for the 
same set of parameters used in Figure|2l but for a frequency 2.0 10^* Hz, in the near-infrared. The model samples lower harmonics and hence 
optically thicker region. Figure [3] shows how diverse can be the curves for different values of h in this regime of optical depths. Therefore 
a set of parameters (as those used in Figure |2j can produce similar curves in optical wavelengths, but not in the infrared. We then consider 
that a 3D approach is more appropriate to model the physical and geometrical properties of the post-shock regions of a polar, particularly 
for mildly optically thick emission. This situation is found in low harmonic numbers or in regions with high densities (A^e > 10^^ cm"''). 
Figure|4]shows the Stokes imaging of the previous models (h = 0.03) in optical and near-infrared frequencies. It illustrates how the total flux 
and polarisation vary along the emitting region. For instance, at phase (0/12) of the lower harmonic, the brighter region (optically thicker) 
is less polarized than the darker region. 

As shock profiles are denser near the white-dwarf surface, it is interesting to check if the above result is also valid for a radially 
variable density and temperature. Figure [5] shows our results for flux and circular polarization in optical and near-infrared wavelengths for 
models having (1) constant or (2) exponentially decreasing (increasing) densities (temperatures). There are differences in the curves which 
are striking in the near-infrared. It is noticeable the raise of the asymmetry in the flux which allow to distinguish the "shock" model from 
a shorter region (compare Figures |5] and [3]l. Also the circular polarization is much larger than when the matter is concentrated near the 
white-dwarf surface. 



^ http://www.hao.ucar.edu/Public/models/pikaia/pikaia.html 
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Table 1. Best fitting to V834 Cen (see text for details). 



Parameter Value 

44 
40 

Ai? 0.132 

AZong (°) 47 

h 0.19 

Si 0.254 

Rth {Rwd) 6.1 

Limits of the threading region in radius (Rwd) 5.32 - 6.9 

Bp (MG) 34.4 

Range of B in the emitting region (MG): 18.8 - 33.2 

Blat (°) 46 

Blong (°) 327 

JVe, maximum electronic density (cm~^) 8.1 10^^ 

T, maximum electronic temperature (keV) and average temperature 171.3 / 20 

Unpolarised component (Jy) 0.0022 (V) - 0.0038 (J) - 0.0026 (H) 

Aphase -0.069 

Attenuation by the upper column 



4 FITTING A DATA SET: APPLICATION TO V834 CEN 



To explore our model and test the fitting procedure of CYCLOPS, we have applied it to multiwavelength observations of a polar (V834 Cen) 

found in the literature. This aims to illustrate the potential of our tool and also to trace the pat h to future improvements. 

V834 Centauri is an object already recognized as having extended emitting regions (e.g.. lCropDej|l989l :l iPotter et al ] |2004l) . The reader 
is referred to these works (and re ferences therein) to a description of the system. From the many published Stokes curves of V834 Cen, 
we have chosen those presented in lBailevetal.1 jl983l) . They represent the magnitudes and circular polarization of V834 Cen in optical and 
near-infrared bands and are ideal to test our fitting procedure using two data set in different harmonic regions. Specifically, we have used the 
V and J band curves converted to flux and polarized flux. We have used frequencies corresponding to the central wavelength of the bands. 

We could not find a good simultaneous fit to both bands using constant values of density and temperature in tangential direction. So we 
adopted a lateral variation as explained in Section l2ni Figure |6] presents the model with the better value of . It was found from a search 
with no restriction in the parameter space. It produces a huge number of possibilities, but we think it is a better approach for testing the search 
algorithm. The parameter of this model is presented in Table [T] The gross features are reproduced, including the relative fluxes between the 
bands. We suggest that a different matter distribution in the emitting region may produce a finer agreement with the data. The above model 
also reproduces the H data of V834 Cen I Bailey et al.lll983 ). even if not included in the fit procedure (Figure|7](. 

The parameters relative to the geometry of V834 Cen a re consistent with previous fittings of the system. For instance, the pre vious 
inclinations found in the literature are in t he range 4 - 6 "ICroppeil I989I: Schwope^et al.j 1993 : Ferrario & Wickramasinghe I990l) . Ex- 
tended emitting region was previous found ('CroppeiJ ll98"i J iFerrario & Wickra masinghe 1990; Pott er et al .Il2004h . This is naturally explained 
by an extended threading region. Pre vious estimates to its radius can be found in Mauche. (200Z , from 5.9 to 37.5 Rwd), using ultraviolet 
photometry and spectroscopy, and in ICroppeil fliiil) . who also found values of order of 30 Rwd using po larimetric data. We can find a 
number of estimates to the temperature using high-energy data: they concentrate in the range 10 - 20 keV dBarlow et al.ll20 06. e.g.). Our 
averag e temperature is also consistent to those values. The magnetic field estimate from Zeeman line splitting in white dwarf photosphere is 



23 MG Ferrario et al. 



1993) . 



The number of possible models, even for narrow bands of all free parameters, is huge. Thus, the search for the best model is limited by 
the available time and processing power. The genetic algorithm has proven to be a robust and efficient method. Even so, a very large database 
of models was generated. It is hard to quantify how much has been explored from the universe of possibilities, but the solution we found is 
an improvement in the sense it simultaneously explain three wavebands. 



5 CONCLUSIONS 

We present a new code, CYCLOPS, to study the optical emission from AM Herculis systems. It calculates the cyclotron emission from a 
3D region on or near the white-dwarf surface. The emitting region is presented by a 3D grid with physical properties spatially variable but 
physically connected by a function. The adopted solution by steps for the radiative transfer is adequate to represent optically thick columns. 
The code also incorporates a procedure to fit an observational data set. An application to previous published multiband V834 Cen data is 
shown. A model with an extended region with electronic temperature and density varying with height and in the tangential direction is able 
to account for the data. 
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Figure 1. (Left hand panel) A diagram representing the white dwarf, the accreting column, and the emitting region (the small white region near the white-dwarf 
surface). It is also shown some important definitions used in the modelling. (Right hand panel) A zoom of the emitting region. The black dot represents the 
anchor point. 
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Figure 2. Examples of Stokes curves for an emitting region with constant electronic temperature and density. The parameters' model are: i = 50°; /3 = 23°; 
B = 29.5 MG; T = 10 keV; Biat = 70°; Bio„g =29°; b= 1.0; AR = 0.20; Aia„a = 90 °. (Left) The soUd (das hed) line was obtained using h = 0.03 (0,10) 
and Ne = 5.13 10^^ (1.55 10^^) cm This model is equivalent to that of lFerrario & Wickramasinghd il99(H . figures 2 to 5, c. We have used a frequency 
of 6.6 10^"' Hz. (Right) The maximum and minimum values of the optical depth among the line-of-sights that compose a image as a function of the orbital 
phase. 
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Figure 3. The same of Figure|2]but using a frequency of 2.0 10^** Hz, in the near-infrared region. 
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Figure 4. Stokes imaging of the models presented in Figures |2] and |3] for h = 0.03. The three lines correspond to the flux, circular polarisation, and linear 
polarisation. The numbers correspond to the system phase. Each line has a differente lookup table where white stands for the highest value. 
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Figure 5. Comparison between models with constant electronic density and temperature and those using "shock-like" profiles. The constant models are 
represented by the solid and long-dashed lines using frequencies of 6.6 10^^ and 2.0 10^'* Hz, respectively. The shock models are the short-dashed and the 
dot-dash lines for the same frequencies. The parameters are those used in the previous figures for h = 0.10. 
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Figure 6. Best fit to the data of V834 Cen from lBailev etalj il983h . The V band is on the left panel and the J band, on the right. See Table[T]for the model 
parameters. 
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Figure 7. Resulting curve in H band using tlie fitted parameters of V834 Cen superposed to the data from lBailev et alj jl983h . 
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APPENDIX A: A NEW EXPRESSION FOR THE SOLUTION OF CYCLOTRON RADIATIVE TRANSFER 



The general solution of the transfer equation for the Stokes parameters in a magnetoactive plasma is given bv lPacholczvk & Swiharll ( Il975[) . 
U nfortunately, there are misprin ts in equations 7 and 13. Therefore, we consider useful to show them again. The correct form of equation 7 
of iPacholczvk & Swiharll ( fl97a> is: 



i = l 
1 

Qh = 'yai{vai+2 + h)Kiey:p{a^s), 

4 

Uh = - ^ exp(ais), 
1=1 

1 

Vh = - - f)Ki exp(Qis). 



(Al) 



The constants Ki i Pacholczvk & Swiharll 1975 . eq. 13) should be written as 

P 



4[(m-n)2 + r2]i/2 
K3A = ~K^-2- 



ai{Io 



I,) - ^^^(Qo - Op) + (t/o - t/p) - ^^^^(v; - V,) 



(A2) 



We refer the reader to the original paper of lPacholczvk & Swiharll ( 119751) for the definitions of the symbols u sed. 



The above expressions were checked as a solution of the differential equation. They also provide the solution of lMeggitt & Wickramasinghe 



19820 for the case with no input radiation (lo = Qo — Uo ~ Vo = 0). 



Now, we present an alternative representation for the complete solution of lPacholczvk & Swiharll ( Il975f) expressed in terms of trigono- 
metric and hyperbolic functions. 



lis) 

Q{s) 

Uis) 

V{s) 
where 

Ki 
K? 



= -2 



2 [K^- cosh (As) + A'f sinh (As) + cos (^ts) + Ki sin (^s)] e""" + 7p, 

I ^ [Kt sinh(As) + iff cosh(As)] + ^ (va^ + h) [-K^ sm(ps) + Kf cos(/is)] | e""" + Qp, 



-2 [a-i [Ki cosh (As) + A'f sinh (As)] + oi [K^ cos (fis) + A'f sin (/is)] } e ''^ + Up, 

2 i^-!l±pi. [K^ sinh(As) + A'f cosh(As)] + ^ (qai - f) [-K^ sin(/is) + cos(/is)] | e""" + 



P 

m 
p 



1 



(lo - h) + [Uo - Up) 



ax 
qa.1 - f 



A 



{Qo-Qp) + ^^^^{Vo-Vp) 



Ki = -JL^tl 



p_fi 
4i? r 



(A3) 



(A4) 



In the above expressions, we have used the following identities: 

Oil = A, 

R = [(m-n)2 + r2]i/2, 

where A, /i, and R are definitions presented in Meggitt & Wickramasinghj Jl982 ). 

In the thermal regime, the particular solution is the Planck function: f Jp, Qp, Up, Vp] = \Bm,0, 0, 0]. 

Equations I A3 1 and I A4I are easily reduced to the solution presented in lMeggitt & Wickramasinghe 1982 ) for the case with no incident 
radiation, [h, Qo, Uo, Vo] — 0. 



